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ABSTRACT 

We study the luminosity function (LF), the comoving rate and the detection rate of Long 
Gamma-Ray Burst (LGRBs) to high redshift, using galaxy catalogues constructed by combin- 
ing high-resolution N-body simulations with semi-analytic models of galaxy formation. We 
assume the collapsar model and different metallicity thresholds, and conclude that LGRBs are 
not good tracers of the star formation history in the universe. Then using the log N — log P 
diagram for BATSE bursts, we determine the LF (with and without evolution with redshift) 
and the formation rate of LGRBs, obtaining constraints on the slope of the power-law. We 
check the resulting redshift distribution with Swift data updated to 2009 August, finding that 
models where LGRBs have as progenitors stars with Z < 0.3^0 and without evolution of 
the LF are in agreement with the data. We also predict that there are about ^ 1% of GRBs at 
redshift z > 6. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are the most luminous explosion in the 
Universe (Zhang & Meszaros 2004). They offer exciting possibili- 
ties for studying astrophysics in extreme conditions, e.g., radiative 
processes in highly relativistic ejecta (Huang et al. 2000; Fan & 
Piran 2008, and references therein). Because of their very large lu- 
minosity, GRBs represent cosmological events, which have been 
detected up to z ~ 8.2 (Tanvir et al. 2009; Salvaterra et al. 2009). 

The observed distribution of the duration of GRBs is bimodal 
(Kouveliotou et al. 1993); long GRBs (hereafter LGRBs) and short 
GRBs, depending on whether their durations are longer or shorter 
than 2 seconds. The current favorite hypotheses for their origin are 
that short GRBs are produced by the merger of compact objects (Li 
& Paczynski 1998; O'Shaughnessy et al. 2008), while LGRBs orig- 
inate from the death of massive stars (with low metallicity), such as 
Wolf-Rayet stars (the collapsar model) (Yoon et al. 2006, 2008; 
Woosley & Heger 2006). Throughout this paper we will deal only 
with the LGRBs. Observational data are consistent with the hy- 
pothesis of the LGRB- supernova connection; at least some LGRBs 
are associated with core-collapse supernovae (Galama et al. 1998; 
Hjorth et al. 2003; Pian et al. 2006; Li 2006; Woosley & Heger 
2006, and references therein). In addition, all supernovae associ- 
ated with GRBs are Type Ic, which supports the hypothesis of 
Wolf-Rayet stars as progenitors of LGRBs. Because of their con- 
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nection with supernovae, LGRBs are potential tracers of the cos- 
mic star formation history (SFH) (Bromm & Loeb 2002; Fynbo 
et al. 2006; Price et al. 2006; Prochaska et al. 2007; Savaglio 2006; 
Totani et al. 2006; Li 2008), nonetheless this connection could be 
non-trivial (see e.g. Kocevski et al. (2009)). 

The detection of LGRBs at high redshift and the connection 
with Type Ic SNe, make them promising for probing the Universe. 
They are probably the only objects that allow us to study the cos- 
mos at high redshift and the early evolution of PopIII stars. There- 
fore, understanding where GRBs are distributed in the universe 
and how they are connected with the star formation rate (SFR) is 
very important. Since the launch of the Swift satellite (Gehrels et 
al. 2004), the number of GRBs with measured redshift has been 
greatly increased. Nonetheless to date there are only ~ 140 Swift 
GRBs with known redshifts. This sample is still too small to con- 
strain their luminosity function (hereafter LF). 

Previous studies (Porciani & Madau 2001; Guetta et al. 2005; 
Natarajan et al. 2005; Daigne et al. 2006; Salvaterra & Chincar- 
ini 2007, and others) touched the problem of determining the red- 
shift distribution of GRBs differently. Their method is based on 
constraining the GRBs distribution by assuming an average energy 
spectrum for all the bursts and that GRBs trace the evolving SFR, 
either with a constant or evolving LF. The redshift distribution, to- 
gether with the LF, can provide important insights not only into 
the physics of the individual objects themselves, but also into the 
evolution of matter in the Universe. 

One of the first pioneering works was Porciani & Madau 
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(2001). The rate of GRB is fitted to observational data, using the 
assumption that SFR is proportional to the GRB rate (Hopkins 
& Beacom 2006). Finally the rate is convolved with a selection 
function dependent on the instrument used. They find that the 
rate of bursts is of about 1-2 GRBs for every one million Type II 
SNe. Notable recent attempts include: Guetta, Piran & Waxman 
(2005) explore a variety of different star-formation rate histories 
and GRB luminosity functions; Natarajan et al. (2005) additionally 
incorporate a simple prescription for a low-metallicity preference 
of GRBs. Daigne et al. 2006 used Monte Carlo simulations to 
predict the GRB evolution, assuming that GRBs follow the SFR, 
the LF is a power-law (independent of redshift) and the peak 
energy is determinated by two relations. They find that: the slope 
of the LF is between 1.5 — 1.7, the Amati relation should be an 
intrinsic relation and, the GRB rate density at 2 = 7 is about 6-7 
times larger than at 2: = 2. They also deduced that the properties 
of GRBs and GRB-progenitors are redshift dependent, since the 
redshift distribution of Swift burst strongly favors their SFR3 
model (see Daigne et al 2006), although that is an unrealistic 
model. 

Salvaterra follows this approach in two different works (Salvaterra 
& Chincarini 2007, Salvaterra et al. 2008), assuming that GRB 
luminosity evolves with redshift and that GRBs form preferentially 
in low-metallicity environments. They use this constraint to set a 
robust upper limit on the slope of bright-end of GRB LF, finding 
that the number of bright GRBs detected by Swift implies that it 
cannot be very steep {5 < 2.6 for progenitors with Z < O.SZq). 
Moreover they found that assuming a threshold of F > 0.4 
[ph/cm^/s], at least ~ 5 - 10% of all detected GRB should lie at 
redshift z > 5. 

In this work, we do not use a GRB comoving rate propor- 
tional to the star formation in the Universe, but only we assume 
that the global rate of GRBs per SNe is on average (over all cos- 
mic times) of about 1 GRB event every 1000 SNe (Langer & Nor- 
man 2006). We derive the LF and formation rate of GRBs using 
a catalogue of galaxies constructed by combining high-resolution 
N-body simulations with a semi-analytic model of galaxy forma- 
tion (Wang et al. 2008). We fit the observed logN-logP relation 
(Kommers et al. 2000) derived from the GRB data of the Burst And 
Transient Source Experiment on board the CGRO satellite (BATSE, 
Fishman et al. 1989) in order to constrain the free parameters of the 
LF. We adopt this method for three GRB progenitor subsamples 
with different cuts in metallicity, following the coUapsar model, 
and assuming a constant and evolving LF. 

By comparing the cumulative distribution of peak photon fluxes, 
Dai (2009) recently proved that the Swift and BATSE samples track 
the same parent population of bursts. For this reason, following Por- 
ciani &: Madau, we rely on the GRBs observed by BATSE as the 
two samples should have comparable LF. 

The paper is organized as follows. We present in section 2 the 
simulated GRB sample used in this work. In section 3, we describe 
the method to reproduce the redshift distribution of LGRB. We de- 
scribe our results in section 4. We give our conclusions in section 
5. 



isi et al. (2009), but we refer to Wang et al. (2008) for a discussion 
of the model. The simulation corresponds to a box of 125 /i^^Mpc 
comoving length and a particle mass of 7.8 x 10* M©. The soften- 
ing length is 5 h^^ kpc. Simulation data were stored in 64 outputs, 
which are approximately logarithmically spaced in time between 
z — 20 and 2 = 1, and linearly spaced in time for 2 < 1. Each 
simulation output was analyzed with the post-processing software 
originally developed for the Millennium Simulation (Springel et al. 
2005). 

In order to extract from the catalog the rate of GRB events, we 
adopt the coUapsar model for LGRBs: all young stars with mass 
> 20 Mq (Larsson et al. 2007) ending their life with a supernova 
should be able to create a BH remnant^. If the coUapsar has high 
angular momentum, the formation of the BH can be accompanied 
by a GRB event (Yoon et al. 2006, 2008). As mentioned in Sec. 1, 
recent studies on the final evolutionary stages of massive stars 
suggested that a Wolf-Rayet (WR) star can produce a LGRB if its 
mass loss rate is small. This is possible only if the metallicity of the 
star is very low. When metallicities are lower than ~ 0.1 — 0.3 Zq, 
the specific angular momentum of the progenitor allows the loss of 
the hydrogen envelope while preserving the helium core (Woosley 
& Heger 2006; Fryer, Woosley & Hartmann 1999). The loss of the 
envelope reduces the material that the jet needs to cross in order 
to escape, while the helium core should be massive enough to 
collapse and power a GRB. 

In order to count the number of GRB events in each snap- 
shot, we select from our catalog objects with redshift between 
< 2 < 9.2, using a procedure similar to that described in section 
3 of Campisi et al. 2009, we count all the possible GRB events in 
the simulated catalog in 3 different subsamples: 
-GRBl, obtained by selecting stars with age < tc = 5 x lO^yr and 
M > 20 M0; 
-GRBl, including stars of age < tc, M > 20 Mq and metallicity 

Zs; 0.3^0; 

-GRB3, defined by selecting stars with age < tc, M > 20 Mq and 
metallicity Z < 0.1^©. 

We compute the number of stars ending their lives as LGRBs, as- 
suming a Salpeter Initial Mass Function (IMF) and that the rate of 
GRB per SNe is on average (over all cosmic times) of about 1 GRB 
event every 1000 SNe (Porciani & Madau 2001; Langer & Norman 
2006)^ 



2.1 Star Formation History 

The coUapsar model links LGRBs to the evolution of massive stars 
whose lifetimes are negligible on cosmological scales. If no other 
condition is required for producing a LGRB event, then the rate of 
LGRBs should be an unbiased tracer of the global star formation in 
the Universe (e.g. Totani 1997; Wijers et al. 1998; Mao & Mo 1998; 
Porciani & Madau 2001; Bromm & Loeb 2002; Fynbo et al. 2006; 
Price et al. 2006; Savaglio 2006; Totani et al. 2006; Prochaska et al. 
2007; Li 2008, and references therein). However, both observations 
and theoretical studies indicate that the metallicity of the progenitor 
star plays an important role in setting the necessary conditions for 



2 SIMULATED LGRB RATE 

In this study, we use the galaxy catalog constructed by Wang et al. 
(2008) for simulations with cosmological parameters from the 
third-year WMAP results. The same catalog was used in Camp- 



1 We also test the case with M > 30 Mq . We obtain that the rate of GRB 
in every box is very close to the rate whit M > 20Mq, since using a 
Salpeter IMF the difference between the two cases is only 0.939 X 10~^. 
^ We test however in Appendix A a different rate (1 GRB every 10.000 
SNe), in order to check how the results change with different assumption. 
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Figure 1. Star Formation Rate History: Log SFR [Mq yr~l Mpc~^] as 
a function of redshift. Dotted line are tlie results from GRB 1 which also 
identify all galaxies in the simulation, dot-dashed line is GRB2 and GRB3 is 
the dashed line. Symbols with error bars are a compilation of observational 
data (Hawkins 2006), the solid line is a best fit of observational data taken 
from Li 2008b. 



a LGRB explosion. In this case, the rate of LGRBs is expected to 
be a biased tracer of the cosmic star formation rate. 

In Fig. 1 we show the SFH for the simulated and for observed 
samples . We compare the cosmic star formation rate obtained us- 
ing all galaxies in the simulation box to that obtained for the three 
samples defined in the Sec. 2 and to the observed SFR in the Uni- 
verse. Since normally the SFR calibrations used for deriving the 
SFH estimates (Hopkins & Beacom 2006) are defined by assum- 
ing the Salpeter (1955) Initial Mass Function, to ensure consistent 
assumptions throughout, we convert SFH estimates in our simula- 
tion'^ to the Salpeter IMF using a simple scale factor. This scale 
factor is established by using the Starburst99 code, which models a 
population-synthesis. 

In Fig. 1, the sample with no threshold on metallicity (GRBl) 
traces the global star formation rate (dotted line in Fig. 1), this 
sample also identify all galaxies in the simulation. This is not 
the case for the two samples with metallicity thresholds (GRB2 - 
dot-dashed line and GRB3 - dashed line). In the GRB2 and GRB3 
samples, the LGRB rate peaks at higher redshift than the cosmic 
star formation rate, as a consequence of the global decrease of 
metallicity with increasing redshift. Data and the best fit (solid 
line) are taken from Hapkins (2006) and Li (2008b). The SFH 
of simulation is in agreement with data within their errors bars. 
As described in Campisi et al. 2009, the peaks of the SFH for 
the GRBs subsamples are shifted at higher redshift, due to the 
selection methods since the objects at higher redshift have lower 
metallicity. 

The important goal of Fig. 1 is to shed some light on the fol- 
lowing issue: do LGRBs trace the SFR? Our results lead to the 
conclusion that LGRBs are not good tracers, albeit they might be 
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biased tracers of SFR. 

But all previous works on the study of redshift distribution of 
LGRBs (see sect. 1) adopted the assumption, i.e. that the GRB 
rate is proportional to the SFR, sometimes convolved with a func- 
tion constraining a metallicity cut-off (e.g Salvaterra & Chincarini 
(2007); Daigne et al. (2006)). We show that selecting in our sim- 
ulation burst with different progenitor's metallicity, the GRB rate 
follows a different evolution with respect to the SFH. Here, the re- 
sults of the work which used a strong correlation between the SFR 
and the LGRB rate could be wrong and should be used with cau- 
tion. 

The results shown in Fig. 1 are in qualitative agreement with 
recent observational estimates (Kistler et al. 2008), and with recent 
theoretical studies also based on the collapsar model (Yoon, Langer 
& Norman 2006 and others). 



3 OBSERVED DISTRIBUTION OF LGRBS 

To predict the observed distribution of the redshift of LGRBs, we 
should take into account that only brightest and pointing toward us 
burst will be observed, so we need to include two important effects: 

I) the coUimation and beaming effects; 

II) the fraction of GRBs seen by the detector (or luminosity func- 
tion of GRBs). 

The number of observed LGRBs is given by: 



A'obs ~ A''rcal fh 



^{L)dL- 



(1) 



where N-^^ai is the total rate of LGRBs in the simulation, /(, is the 
fraction of LGRBs pointing toward us and the integral gives the 
fraction of LGRBs with luminosity bigger than the corresponding 
limit flux of the detector. 



3.1 Beaming effect 

There is a general consensus that GRBs are jetted sources (Waxman 
et al. 1998; Rhoads 1997). This implies fundamental corrections to 
the energy budget and the GRB rates. A canonical GRB does not 
light up the full celestial sphere but rather a fraction, the so called 
beaming fraction (Sari et al. 1998): 

/ft = (1 -cos 61) ^e'^/2 

where 6 is the opening angle of the jet. Thus, the overall GRB rate 
clearly depends on the fact that GRBs are beamed and the rates 
have to be corrected by a factor fb. Thus the true rate integrated 
over a time interval is A^roai = (/t)^^ A'obs. This has been com- 
puted traditionally in terms of the beaming correction factor, which 
is defined as the ratio of the total number of bursts to the observed 
ones. To estimate the overall GRB rate we use the average beaming 
correction /b. The average value of ^ is ~ Qdeg (Ghirlanda et al. 
2007), giving < fb >~ 0.0055. We will use this average beaming 
factor throughout our work. 

3.2 LGRB luminosity function 

The luminosity function (LF) of LGRBs is still poorly tested as 
the data are too sparse for an empirical determination of the burst 
luminosity function. The standard approach to constrain the GRB 
LF from observations (Porciani & Madau 2001, Daigne et al. 2006, 
Salvaterra et al. 2007, and others) is first to assume a model for 
the LF, for the GRB rate, and for the energy spectrum. Secondly 
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the model parameters are constrained by the observed data. To this 
aim it is customary to fit the GRBs observed by BATSE using the 
differential peak flux distribution, logN — logP diagram (Schmidt 
1999). In particular we fit the observed rate of burst with observed 
peak fluxes F between (Fi, F2), described by the equation: 



dN 



iFi<F< F2) = I 
Jo 



dz 



dV{z) An Rg-rb{z) 

dz Att 1 + z 



X / dL'(t){L')e{F), 

'HFi,z) 



(2) 



where (dV{z)/dz) is the comoving volume element'', AQ is the 
solid angle covered on the sky by the survey, the factor (1 + 2)^' 
accounts the cosmological time dilation, i?GRB (2) is the comoving 
GRB rate density, 0(L') is the GRB luminosity function, and e(F) 
is the detector efficiency as a function of photon flux. 

We fit equation (2) using the rate of GRB described in sec- 
tion 2, and we assume the following models for the LF and for the 
energy spectrum. 

3.2.1 Luminosity function 

To model the number of GRBs at different flux limits, we assume 
that the luminosity function has the form: 



*W^i.(^)e.p(-^ 



(3) 



where L is the isotropic equivalent intrinsic burst luminosity, ^ 
is the asymptotic slope at the faint end, L, is the characteris- 
tic cutoff luminosity, and K is the normalization constant so that 
the integral over the luminosity function equals unity. We take 
K = [Lor{-^ - a)]-'^ (for ^ <-l) (Porciani & Madau 2001). 
For the cutoff luminosity we consider differents scenarios, since L* 
could to increase with redshift follow the equation L* — L(z) = 
La{l + zY. 



3.2.2 Energy spectrum 

We assume the empirical form for the GRB spectrum proposed by 
Bandetal. 1993: 
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We adopt the best fit energy spectral indices (i.e., a ^ —1 
and P ~ —2.25) reported in Preece et al. (2000), and the spectral 
break energy distribution at Sb = 511 keV (Porciani & Madau 
2001). These parameters were found by fitting 5500 different spec- 
tra, which is the most extensive GRB sample with spectral charac- 
teristics to date. 

The photon flux F[ph/cm'^/s] observed at the Earth in the 
energy band £„„„ < E < EnmK, emitted by an isotropically 
radiating source at redshift z can be written by: 



F = 



(1 + z) r(i+^'^" 



S{E) dE 



AtvDI 



(5) 



where Dl is the luminosity distance and S{E) is the rest frame 
energy spectrum. We consider £'mi„=50 keV and i5max=300 keV for 
BATSE, while for the Burst Alert Telescope (BAT, Barthelmy et al. 
(2005)) we use iJmin=15 keV and i5niax=150 keV. It is customary 
to define the isotropic equivalent burst luminosity in the rest frame 
photon energy 30 — 2000 keV by: 



^2000 keV 

L= / ES{E)dE, 



(6) 



dV{z)/dz = 4TTcdl{z)/lH{z){l + z)2] 



Thus, combining equations (4-5-6) we get L{F,z) to use in the 
integration limit of eq. (2). 



3.3 Best fit: results 

In our model we have differents parameters to fit: the characteristic 
cutoff luminosity L,, the slope ^ of the LF and a third parameter S 
which is the evolution of L, with redshift, and we fixed it to values 
between 0-3.5 (Salvaterra et al. 2009). We follow the approach of 
Porciani & Madau (2001), using the observed differential number 
counts of BATSE in the range 50-300 keV from Tab. 2 in Kommers 
et al. (2000). The observed sample include 1998 GRBs with peak 
flux in the range 0.18 — 20.0 [ph/cm^/s], and detector efficiency 
described by e{F) = 0.5[1 + erf (-4.801 + 29868 F)] (Kommers 
et al. 2000). Dai (2009) showed that the choice of the BATSE sam- 
ple with respect to the Swift one is equivalent since the two samples 
represent the same population of bursts. Moreover, Dai shows that 
the distribution of the Swift sample matches that of the BATSE sam- 
ple (when approaching the detection limits) so in the follow we can 
use the same trigger efficiency for Swift triggers. 
The fits to the data are done by minimizing the difference of the 
logarithm between model and observational data.This is like a 
simple x^ minimization, but the points are not independent. We 
tried to fit the data directly, finding that the best fit value gave too 
much weight to the central regions. We fitted in logarithm values 
so that the overall shape of the contours has an increased influence 
on the fit. In all considered cases we always find a clear minimum. 
The best-fit model is shown in Fig.2 and the parameters for all 
the subsamples are listed in Table 1. In Fig.2 we show only a com- 
parison between the observed distribution of LGRBs and the pre- 
dicted distribution obtained using the GRB2's subsample, similar 
lines are obtained with the others subsamples and models. The ob- 
served distribution is taken from Kommers et al. (2000), for BATSE 
detector (energy range 50-300 keV). The data are converted into 
rates per unit time per unit solid angle following Kommers's work. 
The horizontal and vertical error bars on data represent the size of 
the energy bin and Poisson uncertainties, respectively. 
Table 1 shows the best fit parameters for all subsamples. The errors 
we quoted are the rms (root mean square) spread in errors from fit- 
ting boostrap catalogue in the same way. 

The given error bars confirm that the characteristic luminosity cut- 
off remains better determined when there is an evolution of the LF 
with redshift. We note that for all subsamples and with very differ- 
ent value of 5, the slope ^ of the LF is well constrained in a range 
of values between 1.686 < — ^ < 1.838 (except one case). These 
values are in agreement with the results in the literature on the slope 
oftheLF, (e.g. in Daigneetal. 2006, where 1.52 < ^ < 1.7). Con- 
versely, the characteristic luminosities span a large range of values 
(0.038 < Lo < 77.875) x 10''°[erg/s], adopting higher values 
when we consider the GRB3 subsample, in particular when the LF 
is constant with the redshift. 
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Figure 2. Comparison between BATSE and best fit model for the log N — 
log P distribution. The dots are the observed BATSE LGRBs in the 50- 
300 keV band (Kommers et al. 2000), and the dark line is the predicted 
distribution with the best fit parameters (we show only for one subsamples, 
similar fits are obtained for the others cases). 
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1.5 
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Sample GRB3 



Lo/(105" ergs-i) 






77.875± 19.243 


2.135±0.184 


1.5 


4.574±0.410 


1.820±0.047 


2.0 


1.970±0.157 


1.787±0.042 


2.5 


0.865±0.064 


1.768±0.039 


3.0 


0.386±0.028 


1.760±0.038 


3.5 


0.173±0.022 


1.760±0.041 



Table 1. Best-fit parameters for different GRBs' subsamples (see section 2) 
and different luminosity evolution with L* = L(z) = Lo(l + z) . Errors 
are computed using boostrap technique. 
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Table 2. Mean redshift for every subsample of GRB with different LP evo- 
lution, compared with the mean Swift redshift. 



4 LGRB REDSHIFT DISTRIBUTION: RESULTS 

About 458 GRBs have been detected by the Swift satellite since 
its launch in 2004 November until August 2009. Among these 
~150 have spectroscopic or photometric redshift determination. 
The number of GRBs with redshift is tightly linked with observ- 
ing conditions, as explained in Jakobsson et al. (2006) (hereafter 
J06). J06 suggested that in order to study the redshift distribution 
of GRBs we should use a subset of all GRBs well placed for optical 
observations. This can be achieved by following 6 criteria, neces- 
sary to "clean" the sample: 1) the burst should have an X-ray posi- 
tion made public within 12 hours; 2) the Galactic foreground should 
be low, i.e. Av < 0.5; 3) the burst should be > 55° from the Sun; 
4) the burst should be not at a polar declination, \dec\ < 70°; 5) 
the burst has to be localised with the XRT; and 6) no nearby bright 
star. Imposing these restrictions does not bias the sample towards 
optically bright afterglows; instead each GRB in the sample has 
favourable observing conditions, i.e. useful follow-up observations 
are likely to be secured. 

Our best-fitting parameters for the LF of GRBs (shown in ta- 
ble 1) are used to predict the redshift distribution for the Swift case. 
We use equation 2 to compute the model prediction of the number 
of GRB with L > L(Liim, z), considering values of F correspond- 
ing to Fiim ~ 0.2 [ph/cm^/s] (in order to compare with Daigne et 
al. 2006). 

In Fig. 3 we compare our model predictions with the number of 
burst detected by the Swift satellite, following the prescriptions of 
the analysis performed by J06 and by using the updated catalog^ 
until GRB 090812. We assume that the observed sample of GRBs 
with redshift determination is representative of all bursts, within 
the error area (Jakobsson et al. 2009). Fig. 3 shows the cumulative 
redshift distribution of observed and simulated GRBs for the 3 dif- 
ferent subsamples. In the model, the expected redshift distribution 
depends on the assumption made on GRB progenitors but also by 
the evolution of the LF with redshift. The distribution in the case 
of the GRB 1 sample is not so far from the observed one. In partic- 
ular it seems that GRB 1 subsample with high evolution for the LF 
(5 > 2.5) reproduce the data. The observed Swift distribution lies 
also close to the distribution of the GRB 2 sample without evolution 
in i,, more at high redshift. This implies that the properties of the 
GRBs do not change with the redshift, since in our simulation the 
progenitor's characteristic does not evolve. 

Figure 3 shows also that the GRB3 subsample is not a good repro- 
duction of observed data, in particular at low redshift (z < 1.5) 
where the lines are outside the error area. 

For completeness in Tab. 2 we also show the comparison be- 
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Figure 3. Redshift distribution for LGRB. Thick step is the observed distribution of Swift burst with sure measured redshift (J06). The blue area takes in 
consideration the error region for the steps, following the procedure of Jakobsson et al. 2009. The upper envelope is produced placing GRBs without redshift 
and those with redshift upper limits at 2: = 0. The lower envelope placing the GRBs without firm redshift at the maximum redshift they can have (giving their 
bluest photometric detection). The model for the expectation of the redshift distribution from our simulation are the dashed lines. For progenitor stars without 
cut in metallicity (GRBl) and with metallicity lower than 0.3 Zq (GRB2) and 0.1 Zq (GRB3). Results are shown for the model with luminosity evolution 
between 0-3.5. 



tween the average value of the redshift for the updated GRBs' sam- 
ple of J06 and the simulated one. We predict that the value of the 

< z > evolves with the threshold in metallicity, in agreement 
with the SFR evolution for the three subsamples. In fact, while the 
GRB 1 prediction gives an average redshift between 1.71 — 2.33 (for 
different values of 5), GRB2 and GRB3 have higher values up to 

< 2 >~ 3.44 for the extreme case where 5 = 3.5 and^ < O.IZq. 
However the observed LGRB have an average redshift of < z >~ 
2.27, which is more lower than the one predicted from the GRB3 
model. From Fig.3 and Tab.2, we are able to rule out the GRB3 
subsample and we conclude that either GRBl with evolving LF 
and GRB2 with a non-evolving LF are possible model within error 
bars. 



4.1 Bright and Faint LGRBs 

Also if there would be an evolution of GRB LF with time (high- 
z GRB are typically brighter than low-z ones) from Fig.3. seems 
to need a very strong evolution to reconcile GRB 1 with observed 
data (S > 2.5). This evolution should imply that the properties of 
GRBs change more with the redshift. This scenario seems to be 
unrealistic. 

In order to step over the detection problem, following the 
same approach of Daigne et al. (2006), we defined two subsam- 
ples of Swift bursts by selecting those with peak flux Fwm > 1-0 
[ph/cm^/s] (Bright) and those with i^i,„, < 1.0 [ph/cm^/s] (Faint). 
In Fig.4 we show the cumulative redshift distribution of observed 
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Figure 4. Redshift distribution for Swift observed and expected burst for 
Faint burst with peak flux F ^ 1 [ph/cm^/s] and BrigAtLGRB with F > 1 
[pli/cm^/s]. The lines are the simulations and the steps are the Swift data. 



Rate 


Faint 


Bright 


% 


z > 5 2 > 6 


z > 5 z > 6 


GRB2 


-3.8 ~ 1.1 


~ 2.2 ~0.9 


Obs 


-8.5 ~2.8 


~4.4 ~2.9 



Table 3. Fraction of GRBs with redshift z > 5 for GRB2 subsample with- 
out evolution for the LP, and observation data from J06, for burst with peak 
flux F ^ 1 [ph/cm^/s] and with F > 1 [ph/cm^/s]. 



and expected bursts. We decided to show only the GRB2 subsam- 
ple, without evolution in LF, since it seems to be the more realistic 
model from Fig.3. 

It is evident that Bright and Faint objects have different redshift 
distribution, both for expected and observed ones. We note that 
the Brigiit observed sample almost overlaps the simulated distribu- 
tion of bursts, overcoming the problem of Fig.3. Conversely, there 
are more Faint bursts at low redshift than predicted in the simula- 
tion. Both Brigiit and Faint have similar behavior a redshift z > &. 
However we argue from Fig.4 that GRB2 subsample with a non- 
evolving LF is the best possible model to reproduce Swift data. 

We quantify the probability to find burst a redshift 2 > 5 and 
a > 6 in Tab. 3 for both subsamples. We expected from model to 
have about 1% of GRBs with redshift z > &. For both cases we find 
results lower than Daigne et al. (2006), which found that at z > 6 
the fraction of Bright bursts should be ~ 2 — 6%. The difference 
between our estimation and Daigne's one is that their SFR3 (pre- 
ferreed by redshift distribution of observational data) is probably 
unrealistic as written in the conclusion of Daigne's paper, suggest- 
ing that their results provide strong evidence that the properties of 
GRBs are redshift dependent. Instead in our case we take into con- 
sideration of the GRB luminosity function by considering a model 
where the characteristic luminosity does not change with redshift, 
with the progenitor properties unchanged. 

We claim that if the distribution of the observed sample will 
not change with the increase of the number of observed LGRB, the 
LF should not have an evolving cutoff among our modelled values. 



Moreover the correlation between Swift bursts and GRB2 sample 
implies that LGRB should be produced by a very massive star with 
metallicity Z < 0.3^©. 



5 DISCUSSION AND CONCLUSIONS 

In this work, we have studied the luminosity function (LF), the 
comoving rate and the detection rate of Long Gamma-Ray Burst 
(LGRBs) in the context of a hierarchical model of galaxy forma- 
tion. 

Assuming the collapsar model and imposing different metal- 
licity constraints we find that: 

• GRBs with low metallicity progenitors (Z < 0.1 — 0.3Zq) 
do not represent a perfect tracer of the cosmic star formation 
history (see Fig. 1). The deviation of the LGRB rate from the star 
formation rate decreases with increasing redshift (as a consequence 
of the global decrease of metallicity with increasing redshift) and 
the bias is stronger as the metallicity threshold assumed is lowered. 

• The LF of LGRBs is well descripted by a power-law with 
exponential cut-off, with well determined slope between 1.686 < 
(, < 1.838. Conversely, the characteristic luminosity spans a large 
range of values (0.038 < io < 77.875 x 10^''[erg/s]), increasing 
at lower metallicity threshold (see table 1). 

• It is possible to reproduce, within error bars, the redshift 
distribution of a subset'' of Swift data with F > 0.2 [ph/cm^/s], 
using: (I) a model without cut in metallicity with a very strong 
evolving LF, or (II) a model with metallicity threshold Z < 0.3^q 
and a non-evolving LF (see Fig.3). Selecting only Bright LGRBs 
(F > 1 [ph/cm^/s]) and suspecting that a scenario where the 
properties of GRB change so strong with redshift seem to be 
enough unrealistic, we rule out the (I) model (see Fig. 4). 



• We predict to have 
{z > 6). 



1% Bright bursts at high redshift 



Our work constrains the LF function of LGRBs, using a rate 
of GRB not proportional to SFR, and assuming different metallicity 
for the progenitors, giving us the possibility to assert that the most 
probable model for LGRB's progenitor have Z < O.SZq and no 
evolution in the LF. 
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APPENDIX A: GRB RATE 

In order to compute the rate of GRB in our simulation, in section 
3 we assume to have 1 GRB every 1000 SNe globally in the Uni- 
verse, in agreement with Porciani & Madau (2001); Langer & Nor- 
man (2006). Nevertheless we explore also how change the results 
assuming a different rate for GRBs. This is what we want test in 
this appendix. We repeat all the work, using a rate of 1 GRB every 
10000 SNe. We fit the BATSE sample, obtaining of course differ- 
ents values for the luminosity function, and we will use this value to 
compute the redshift distribution. To compare with Eig.3, we show 
in Eig. Al the redshift distribution for GRBl and GRB2 subsam- 
ple, for GRB3 subsample the rate of GRB is too low to find a best 
fit for the LE of the BATSE data. 

The dependence of the LE with the redshift is not remarkable, since 
(as in GRB3 in Eig.3) the rate of GRB is smaller than previous case. 
In Eig. Al it is evident that the best model to reproduce the data is 
the GRBl subsample (with and without LE evolution). We expect 
this results because the assumption to have 1 GRB every 10000 
SNe calls for a big number of SNe where there are also more low 
metallicity objects, and for this reason the redshift distribution of 
GRBl and GRB2 is shifted at high redshift respect Eig.3. 
We conclude that, since models and observations suggest that a 
metallicity dependency is required, a rate smaller than that using 
in our work could be unrealistic. 



© 2010 RAS, MNRAS 000, 1-7 



Redshift distribution and luminosity function ofLGRBs 9 



GRB-1 



GRB-2 



a 
o 



fc 




N 

V 



a 
o 



fc 




2 3 4 5 6 7 
Redshift 

Figure Al. Redshift distribution for LGRB. Thick step is the observed distribution of Swift burst with sure measured redshift (J06). The blue area talces in 
consideration the en'or region for the steps, following the procedure of Jakobsson et al. 2009. The upper envelope is produced placing GRBs without redshift 
and those with redshift upper limits at z=0, instead the lower envelope placing the GRBs without firm redshift at the maximum redshift they can have (giving 
their bluest photometric detection). The model for the expectation of the redshift distribution from our simulation are the dashed lines. For progenitor stars 
without cut in metallicity (GRBl) and with metallicity lower than 0.3^0 (GRB2). Results are shown for the model with luminosity evolution between 0-2.5. 
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